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Abstract 

We present a new formalism for the theory of relativistic dissipative hydrodynamics. Here, we 
look for the minimal structure of such a theory which satisfies the covariance and causality by 
introducing the memory effect in irreversible currents. Our theory has a much simpler structure 
and thus has several advantages for practical purposes compared to the Israel-Stewart theory (IS). 
It can readily be applied to the full three-dimensional hydrodynamical calculations. We apply our 
formalism to the Bjorken model and the results are shown to be analogous to the IS. 
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I. INTRODUCTION 



The ideal hydrodynamical description for the dynamics of hot and dense matter achieved 
in RHIC experiments works amazingly well, particularly for the behavior of collective flow 
parameters. Together with other signals, the success of the approach is considered as the 
indication of the emergence of a new state of strongly interacting matter, the plasma of 
quarks and gluons (QGP). The comparison between RHIC and SPS results shows that this 
new state of matter is formed at the very early stage of the relativistic heavy ion collisions 
for RHIC energies. The nature of the QGP seems rather a strongly interacting fluid (sQGP) 
than a ideal free parton gas [1] which still flows like an ideal-fluid. This idea, that the QGP 
behaves as real ideal fluid, raised an interesting perspective, since the viscosity for the strong 
coupling limit of 4D conformal theory obtained from the supersymmetric string theory in 
an (Anti de Sitter) AdS space found to be very small On the other hand, Hirano and 
Gyulassy argue that this is due to the entropy density of the QGP which is much larger 
than that of the hadronic phase (3]. 

At any rate, we know that there still exist several open problems in the interpretation of 
data in terms of the hydrodynamical model [4| . These questions require careful examination 
to extract quantitative and precise information on the properties of QGP. In particular, we 
should study the effect of dissipative processes on the collective flow variables. Several works 
have been done in this direction [5j. However, strictly speaking, a quantitative and consistent 
analysis of the viscosity within the framework of relativistic hydrodynamics has not yet been 
done completely. This is because the introduction of dissipative phenomena in relativistic 
hydrodynamics casts difficult problems, both conceptual and technical. Initially Eckart, and 
later, Landau-Lifshitz introduced the dissipative effects in relativistic hydrodynamics in a 
covariant manner 

aa 

. It is, however, known that their formalism leads to the problem of 
acausality, that is, a pulse signal propagates with infinite speed. Thus, relativistic covariance 
is not a sufficient condition for a consistent relativistic dissipative dynamics 

To cure this problem, relativistic hydrodynamics in the framework of extended thermo- 

n nnn 

dynamics was developed by Miiller [10] and later by Israel and Stewart [111, [12|, [13[ . From 
the kinetic point of view 14j, this formalism corresponds to the extension of equilibrium 
thermodynamics to include the second order moments of kinetic variables. This is the reason 
why this theory is usually referred to as the second order theory. However, this theory is too 
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general and complex, introducing many unknown parameters from the point of the present 
QCD dynamics. Moreover, the irreversible currents are treated as independent variables in 
addition to the usual hydrodynamical degrees of freedom (the velocity field plus densities of 
extensive quantities) and should be obtained by solving highly coupled partial differential 
equations. These features make the application of the theory to practical problems very 
difficult. 

On the other hand, even if we overcome the technical difficulties to apply the complete 
theory to the analysis of relativistic heavy ion reactions, it is very unlikely that we can 
extract such a detailed information from the experimental data at present moment. This 
is because there exist many uncertainties of the hydrodynamical approach itself such as 
the initial condition, finite size correction, event-by-event fluctuations, particle decoupling 
process etc. |4|. However, there exists some interesting questions which requires a correct 
treatment of the viscosity in a realistic 3D simulations, such as a possible generation of shock 
wave and its propagation in the QGP triggered by a jet[15|- Analysis of such phenomena 
will certainly give important information of new, genuine hydrodynamical properties of the 
QGP. Therefore, what we need now is not a very general theory of relativistic dissipative 
hydrodynamics with so many unknown parameters but one with the minimum number 
of parameters necessary to preserve fundamental principles such as causality and entropy 
production. 

As a matter of fact, it is known that the Israel-Stewart theory (IS) is not the unique 
approach to relativistic dissipative hydrodynamics. To authors' knowledge, there is at least 
one other causal theory called the divergence type 0, 0, Q, Q, [20j. In this work, we 
present a very simple alternative theory which satisfies the minimal conditions mentioned 
above. We show that the causality problem of the Landau-Lifshitz formalism can be solved 
by introducing a memory effect. This memory effect is characterized by the relaxation 
time tr, so that our theory introduces only one additional parameter to the usual bulk 
viscosity, shear viscosity and thermal conduction coefficients of the Navier-Stokes equation. 
Our theory recovers the relativistic Navier-Stokes equation in the limit of vanishing values 
of this relaxation time. 

As described later more in detail, our approach has a fundamental advantage from the 
practical point of view in addition to its physical simplicity. The dissipative terms are 
explicitly given by the integral of the independent variables of the usual ideal hydrodynamics. 
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Thus, the implementation of our method to the existing ideal hydro-codes is straightforward, 
oarticularly, to those based on the local Lagrangian coordinate system such as SPheRIO 
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22j. 



The present paper is organized as follows. In the next section, we briefly review the 
problem of the acausal propagation in the diffusion equation and the method to cure this 
problem in terms of the memory effect, which leads to the so-called telegraphist's equa- 
tion. In the Sec. Ill, we analyse the structure of the Landau-Lifshitz theory of relativistic 
dissipative hydrodynamics and introduce the memory effect to solve the acausal problem 
due to its parabolic nature. We thus obtain the dissipative hydrodynamical equations with 
the minimum number of parameters which satisfies causality. In Sec. IV, we discuss the 
problem of entropy production in our formalism. In Sec.V, w e apply our equation to the 



Bjorken solution, and compare with the previous analyses 
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23| . In Sec. VI, 



we summarize our results and discuss possible immediate applications. 



II. DIFFUSION EQUATION AND ACAUSALITY 

The fundamental problem of the first order theory like the Navier-Stokes theory comes 
from the fact that the diffusion equation is parabolic. This means that the velocity of signal 
propagation is infinite. This problem was first addressed by Cattaneo [28] in the case of 
heat conduction. He discussed that the problem of acausal propagation in usual diffusion 
equation, 

|n = CV 2 n, (1) 
can be cured by the introduction of the second order time derivative as 

d 2 n dn 

thus converting the parabolic equation to a hyperbolic one. Here, ( and tr are a diffusion 
constant and a relaxation time, respectively. Eq.p]) is sometimes called as telegraphist's 
equation. For a suitable choice of the parameters tr and (, we can recover the causal 
propagation of diffusion process. In fact, the maximum velocity of the signal propagation 



governed by this equation is shown to be 
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To be consistent with the special relativity, we should have 

r R > J, (4) 

The diffusion equation, Eq.([T]), corresponds to tr = and hence v max — > oo. 

The physical origin of above telegraphist's equation can be understood as follows. The 
diffusion process is a typical relaxation process for a conserved quantity. Thus, it should 
satisfy the equation of continuity, 

| + V-J=0. (5) 

In non-equilibrium thermodynamics, the irreversible current j is assumed to be proportional 
to a thermodynamic force, 

j = -LF, (6) 

where the Onsager coefficient L is, in general, a function of thermodynamic quantities. Here, 
we assume it to be constant for the sake of illustration. When n is a number density, the 
thermodynamic force is given by the gradient of n because of Fick's law, 

F = Vn. (7) 

Substituting Eq.flHJ) into Eq.flHJ), we get the diffusion equation (TjQ). 

Fick's law tells us that the above diffusion process is induced by inhomogeneous distribu- 
tion. In Eq.flE]), the space inhomogeneity immediately gives rise to the irreversible current. 
However, this is a very idealized situation. In general, the generation of irreversible currents 
has a time delay. Thus, we may think of memory effects within the linear response of the 
system. Phenomenologically, this can be done by introducing the following memory function 



G (tX) = —e-( t - t 'V TR ,t>t' (8) 
= 0, t<t' (9) 

where tr is a constant (relaxation time) and rewrite Eq.([n]) as 

]=- f G(t,t')LF(t')dt'. (10) 

J — oo 
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In the limit of r# — > 0, we have G(t,t') —> 5 (t — f) so that the original equation is 



recovered 



31] . The time derivative of the irreversible current fflUj) leads to 



| = _±^ (i) _±j. 

Tr Tr 

This is called the Maxwell-Cattaneo type equation. Substituting into the equation of conti- 
nuity (j5j), we arrive at 

d 2 n „ / 1 ^ . , 1 - 
1 / dn 



Tr V 

which is exactly Eq. (j2J). 

We here showed the problem of acausality in the diffusion equation, which is based 
on non-equilibrium thermodynamics (the first order theory), and how this can be solved 
by introducing the memory function. For more microscopic derivation of telegraphist's 
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equation and foundation of the second order theory, see Refs. 0, Q, 
it was recently shown that the macroscopic equation of motion, which is obtained by using 
systematic coarse-grainings from the Heisenberg equation of motion, is not the diffusion 
equation but telegraphist's equation j^ . 

Till now, we have emphasized the importance of the memory effect in the view of causality, 
but there exists another important reason to introduce the memory effect. It was pointed 
out that the diffusion equation also contradicts with sum rules associated with conservation 
aws. This can also be solved by replacing the diffusion equation with telegraphist's equation 
32|. Relativistic dissipative hydrodynamics should be consistent with QCD, where there 
are conserved quantities and associated sum rules. In this sense, the presence of memory 
effect should be inevitable. 



III. MINIMAL CAUSAL HYDRODYNAMICS 

Eckart and Landau-Lifshitz derived the relativistic dissipative hydrodynamics following 
non-equilibrium thermodynamics as discussed in the preceding section p, Q]. Their the- 
ories are just the covariant versions of the Navier-Stokes equation and the corresponding 
equations still continue to be parabolic. Therefore, they do not satisfy causality and some 
modification should be required. In the IS and the divergence type theory, the definition 
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of the entropy four-flux is generalized and, to satisfy the second law of thermodynamics, 
modified thermodynamic forces are obtained. In this section, we propose another approach, 
where the problem of acausality is solved by introducing the memory effect like Eq. (11 01) . 

For this purpose, let us first analyse briefly the structure of the Landau-Lifshitz theory 
(LL). The hydrodynamical equation of motion is written as the conservation of the energy- 
momentum tensor , 

dyT^ = 0, (12) 
together with the conservation of a quantity, for example, the baryon number, 

d^N" = 0. (13) 

In the LL, it is assumed that thermodynamic relations are valid in the local rest frame of 
the energy-momentum tensor. The energy-momentum tensor is expressed as 

= euV - {p + U) + 7f' MU , (14) 

where, e, p, IT and are respectively the energy density, pressure, four velocity of the 
fluid and bulk and shear viscous stresses. In the LL, the velocity field is defined in such a 
way that the energy current vanishes in the local rest frame, —>■ (1,0,0,0). In this local 
rest frame, it is assumed that the equation of state and thermodynamical relations are valid 
as if the fluid were in equilibrium. As usual, we write 




where 7 is the Lorentz factor and 

m% = 1. 



The tensor P^ u is the projection operator to the space orthogonal to VtF" and given by 

which has the following trace property, 

f>A A = £- (15) 

A=0 

Here, D is the number of spatial dimensions. 
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In the LL, the current for the conserved quantity (e.g., net baryon number) takes the 
form 

iV M = rax^ + iS, (16) 



where is the 
the constraints 



leat conduction part of the current. For the irreversible currents, we require 

= 0, (17) 



11C 
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and 

uX = 0. (18) 

These constraints permit us to interpret e and n respectively as the energy and net baryon 
number densities in the local rest frame. In fact, from Eq. (1181) . in the rest frame, we have 



so that n is the net baryon number density in the local rest frame. It is noteworthy that, 
different from the ideal case (or Eckart case), iV° is not related with the net baryon number 
density n in any frame by the Lorentz contraction, but 

N° = wy + [p ■ v) 7. 

With these irreversible currents, of course, the entropy is not conserved. Instead, from 
Eqs.(TT2l) and ( fT3l) with the constraints Eqs. ( TT7j) and ( TT8l) . we have [7|] 

(su^ - av") = i (-P^U + tO d„u v - u^a, (19) 

where a = fi/T with temperature T and chemical potential fi. Landau-Lifshitz identifies 
the term 

= su^ - (20) 

as the entropy four-flux. The r.h.s. of Eq. (Tl9|) is the source term for entropy production. 

In non-equilibrium thermodynamics, it is interpreted that entropy production is the sum 
of the products of thermodynamic forces and irreversible currents. Thus, we can define the 
scalar, vector and tensor thermodynamic forces, 

F = d a u a , F^ = <9 M a, F^ = d^, 



S 



respectively. To satisfy the second law of thermodynamics, we assume that the entropy 
production is positive, 

_L (_p^ n + vrH d,u u - u^a > 0. (21) 

To maintain this algebraic positive defmiteness, the most general irreversible currents are 
given by linear combinations of the thermodynamic forces with the coefficients appropriately 
chosen. However, if we accept the Curie (symmetry) principle which forbids the mixture of 



different types of thermodynamic forces 34J |. the irreversible currents are given by 



n = -(F = -(d a u a , 

v„ = P^v v = -kP^F" = -nP^d u a, (22) 

where (, rj and k are bulk viscosity, shear viscosity and thermal conductivity coefficients, 
respectively. Here, P^ a ^ is the double symmetric traceless projection, 

pfJ.va/3 _ }_ ^ptiapvfl _|_ pfi/3 pua^j ^_p^u pa/3 (23) 

and we have introduced the quantities TT al3 and v v which correspond respectively to the shear 
tensor and irreversible current before the projection. They are in general not orthogonal to 

so the projection operators are necessary to satisfy the constraints Eqs.(fT71) and ( 1181) . 

Eqs. (l22|) are the prescription of the LL. As mentioned, the LL leads to the acausal prop- 
agation of signal. So we should modify these equations to satisfy the relativistic causality 
principle. The basic point is that the equations of the LL form a parabolic system and we 
have to convert it to the hyperbolic one. However, at this moment, the generalization of 
these equation in order to obtain hyperbolic equations is rather self-evident. We introduce 
the memory function in each irreversible currents, Eq. (1221) . 



n(r) = - r dr'G(r,T')Cd a u a (r'), 

J — oo 

5r^ (r) = / dr'G (r, r') rjdW (r') , 

dr'G (r, r') nd^a (r') , (24) 
where r = r (r, t) is the local proper time. As before, the shear tensor ir^ u and the irreversible 
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current z/ M are then given by the projection of these integrals as 

tt^ = P^K af3 (r) , 

When we start with the finite initial time, say, To, the above integrals should read 

n (t) = - / dr'G (r, r') (^0 a u a (r') + e- (r - ro)/rB n , (26) 

^ TO 

tt"" (r) = / dr'G (r, r') rjdW (r) + e - (r - To)/Tfl F%, (27) 

J TO 

z7" (r) = - / dr'G (r, r') (r) + e^^^o, (28) 

J TO 

where IIo, 7r^o and z/^o are initial conditions given at tq. The bulk viscosity can be deter- 
mined by T^ v (r ) and (r ) as 

no = ^(^o-T;(r ))-po. 

As for tt^o and Z/^o we observe that they are the projected part of T^ u (tq) and (tq), 
respectively, 

^o = ^ (r )-n < = P^(r ), 

tt^o = (r ) - (e + Po + n ) « + (po + n )<T = P^T^{t ). 
This is equivalent to set up the following initial conditions, 



^0 = ^ (r )-7io<, 

tt^o = (r ) - ( £o + Po + n ) + (po + U )g^. 

The initial four-velocity Uq is determined as the time-like eigenvector of the matrix T^ v (r ) 
with eigenvalue e . Using the equation of state together with n = u^N 11 , we can determine 
the pressure po. 

In Eqs. fl24l) . we have used the same memory function G and consequently a common 
relaxation time tr for the bulk and shear viscosities and heat conduction. We could have 
used different relaxation times for each irreversible current and this would not alter the 
basic structure of our theory. However, here we stay with a common relaxation time for 
all of them for the sake of simplicity. We consider the situation where the time scales of 
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the microscopic degrees of freedom are well separated from those of the macroscopic ones. 
Then, the effect of the differences of the microscopic relaxation times should not be much 
relevant in the dynamics described in the macroscopic time scale. Thus we just represent 
these microscopic time scales in terms of one relaxation time tr. 

The integral expressions (1241) are equivalent to the following differential equations, 

TT ^ « dU 

CLT 

drc^ 

ft" = rjO^u" - t r ^—, 
dr 

W = -nd^a - t r ^ } (29) 

GST 

where 

is the total derivative with respect to the proper time. In a practical implementation of 
our theory, we may solve the above differential equations together with the other hydrody- 
namical part. In this case, if we wish, we may add also the terms which violates the Curie 
principle, without any extra difficulties. The above equations, after the projection (1251) . can 
be compared to the corresponding equations in the simplest version of the IS, which are 
obtained phenomenologically based on extended thermodynamics, 

11/5 = ~\Cis (d a u a + Po-^- - a d a u c 

/ dn^ \ 

n% = 2 VIS P^ (d a u p - /3 2 ^ - a x d a (u IS )A , 

v? s = -kjsP^ (-^p9 v a + Pi^ + aod v U IS + a,d a (n IS )A , (30) 

where a , ati, /3q, (3\ and (3 2 are constants. Note that the definitions of parameters r), ( and 
k are different from that of the IS. Eqs. (1291) and (1301) have similar aspects, in particular, 
if we include the terms which violate the Curie principle and introduce the 3 different 
relaxation times. However, our equation is not a special case of the IS. In the IS, the 
projection operators, which are necessary to satisfy the orthogonality conditions ( fTTl) and 
f fT8|) . are included in the differential equations themselves. Thus, it is not possible to derive 
our equations from the IS. For example, in our theory, we can write down the differential 
equation of the heat conduction by using Eq. ([29]) as follows, 



V* = -KP^dvCL - Tr — + 

dr dr 
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The last term of the above equation do not appear in the IS. 

In practice, the differential equations for the irreversible currents can be solved together 
with the equations of motion of the usual hydrodynamic variables in a coupled way. In the 
case of the IS, due to the presence of projection operators, these irreversible current should 
be determined simultaneously with the equations for acceleration of the fluid. Thus, in the 
general 3D case, (14 x 14) matrix inversion is required in each time step. In our case, due 
to the memory effect integral, the equation for the acceleration of fluid can be determined 
only from the past values of the irreversible currents, so that we need the inversion of at 
most (3 x 3) matrix for the acceleration of the fluid. The time derivatives for the irreversible 
currents are decoupled from the normal hydrodynamic variables. 

In spite of these differences, our equations are found to be still hyperbolic in the linearized 
form. When we consider the propagation of small perturbations on the homogeneous and 
static background, the projection operator turns out to be a constant matrix. Therefore we 
can easily see that our linearized equation of motion has the same structure as the IS with 
an = a>i = 0. Thus the speed of pulse propagation is finite as discussed by Hiscock-Lindblom 



35|,l36|,l37|. 



In our case, we can explicit the space component of the four- divergence of the shear tensor 
TT flu can be written as 



d 1 

dr a 



7 3 ' 



u 



+ 



dn^ l )| 



dr 



1 d M*7r2 



(31) 



where we have introduced the symmetric tensor, 



2 v 



a/3 



+ ^) , 



and used the three-vector notation, 



/ ^(a 1) \ 



71 



7T V 



7T 



(a 2) 



. 3) 

Furthermore, for the sake of later convenience we have introduced the "reference density" 
a = a (t, r) defined by 

1 rlrr 

(32) 



1 da 
a dr 
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This will play an important role for the smoothed particle hydrodynamics (SPH) formulation 



21 



221 ] . Note that, in general for any function / = / (t, r) we have the identity, 



dr \a 

The equation of motion finally has the form 



d 1 
dr a 



e + V + n - ^u a ^ a °) + ^a + l (u a u^ a P) ) } u l - u a ^ a 

= —ft ( P + n + i - i^ifi,T + iv • (tt q - . (34) 

Note that the right-hand side of the above equation does not contain the acceleration. From 
this equation, we can determine the time derivative of u\ i — 1,2, 3. 
The time component is equivalent to the energy (entropy) equation, 

d^su 11 ) = - (-P^n + tt^) <9 M ^ + adX- (35) 

so that using Eq. (l33|) . we have 

^- (-) = != {-P^U + Tin d,u u + ady. (36) 

GET \0 J (TT 

The conservation of net baryon number is written as 

A. (2) = _ V». (37) 
dr \<jJ a 

Eqs. f )32|[34ll36ll37p specify the time evolution of u l , s,n and a. These equations constitute a 
closed system with the use of the equation of state and thermodynamic relations, fi = de / dn 
and T = de/ds, and also the integral expressions, Eq. (j24p for II, tt^ v and z/ M (and their 
projected tensors, tt^ u and v^). 

From the practical point of view, our approach has several advantages to the IS in addition 
to its simple physical structure. Although the IS invokes the information from the kinetic 
theory, it is not probable that we can determine the parameters from the relativistic heavy 
ion processes quantitatively at the present moment. From the theoretical point of view, 
one could try to fix these parameters by using the Boltzmann equation approach to QCD. 
However, the Boltzmann equation is applicable to describe the behavior of a gas and if 
the QGP is a strongly interacting fluid, then such an approach will fail. Rather, we prefer 
to keep the simplest physical structure of the Navier-Stokes theory preserving causality. 
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To accomplish this, we simply introduced the memory effect characterized by the relaxation 
time. All irreversible currents are then expressed by Eg. (1241) . We only need the past values of 
the independent variables of the usual ideal fluid dynamics to advance in time. The integral 
expressions (|24p are easy to be evaluated when we use the Lagrangian coordinate system. 
Because of the simple form of the dissipative terms, it is straightforward to incorporate these 



eguations to the realistic full three-dimensional hydro-code such as SPheRIO [2ll. |22|. 

Until now, we have considered that the relaxation time tr is constant. However in 
practical problems, it is a function of thermodynamical variables. Then the memory function 
should be generalized as 

G(T,rO-^-e- /; '^ r ''. (38) 
Even in this case, Egs. (l24l) and (134"]) are still valid. 

IV. ENTROPY PRODUCTION 

It should be emphasized that our theory is not a simplified version of the IS but there 
exists an essential difference for the treatment of the entropy production term. The IS 
reguires the general algebraic form of the non-negative definite expression for entropy pro- 
duction following non-eguilibrium thermodynamics. In our approach, we have relaxed this 
condition, that is, the expression Eg. ( 1211) for the entropy production 

I (_p/-n + tO d^u v - u^a > 0, 

does not guarantee algebraically owing to the non-locality in time contained in II, ir^ and z/ M 
through Eg. (1241) . This might seem to be dangerous. However, strictly speaking, the increase 
of entropy is essentially a concept in the eguilibrium thermodynamics and the reguirement 
of positiveness should apply only to thermal eguilibrium states. As a matter of fact, it 
was recently shown that the entropy absorption process can occur in the non-eguilibrium 



processes of mesoscopic systems 
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In our approach, we are dealing with a fluid element 
which is out of eguilibrium, interacting with the neighboring elements. Therefore, within the 
relaxation time, its entropy content may increase or decrease depending on the dynamics and 
its time scales. Thus, the reguirement of the algebraic positive definiteness irrespective of 
any field configuration seems to be too restrictive. The reguirement of non-negative entropy 
production may be relaxed for far from eguilibrium states. In our case, apart from the 
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projection operators, the expression for entropy production has the form 

1 



Q (r) = / (r) 



Tr 



J T dr'e-^- T '^ TR f{T r ) 



(39) 



where / is one of c^ttj, or d^a. For small tr, we may expand / (V) near r, 



/ (r') = / (r) - (r - r') 



n df (r) 



+ 



and we have 



Thus, as far as 



Q(r) = f (r) 



f (r) - r R 



df{r) 
dr 



O (r R 2 ) 



df(r) 



dr 



<\f(r) 



(40) 



(41) 



(42) 



the positiveness of the entropy is ensured. The l.h.s. of the above equation is the amount 
of variation of / (r) within a small time interval tr. Thus, the above condition shows that 
if the time evolution of the system is not too violent (the change of field values within the 
relaxation time is less than its value), then the local entropy production is not negative. 
For the example discussed below, we can show explicitly the positive definiteness of entropy 
production on our formulation. 

In the above, we considered the relaxation time as constant just for the illustration. The 
similar conclusion can be derived when the variation of the relaxation time is not so violent. 



V. ONE-DIMENSIONAL SCALING SOLUTION 



To see how the above scheme works, let us apply it to the one dimensional scaling solution 



of the Bjorken model. This has been studied already in the framework of the IS 
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, [27|] . The components of the irreversible currents is written down explicitly as 

II (r) = - f dr'G (r, r') ^ + r R (r ) G (r, r ) n (r ) , (43) 

J TO T 

fi (r) = - f dr'G (r, r') ^ + (r ) G (r, r ) Q (r ) , (44) 
$ (r) = - jf dr'G (r, r') «^ + (r ) G (r, r ) $ (r ) , (45) 

( — sinh 2 y sinh y cosh y \ 
2 n(r), (46) 
sinh ?/ cosh ?/ — cosh y I 

cosh u \ 

U(r), (47) 
— sinh ?/ / 

where we have used the hyperbolic variables, 

t = r cosh y, x = t sinh y, 

and used the scaling ansatz in y (that is, there is no y dependence in thermodynamic 
variables). II (tq) , Q (tq) and $ (to) are initial values for II (r) , Q (r) and $ (r) . We obtain 

9Q 

7T^ = P^Tfap = P"", (48) 

3 

= = 0. (49) 

As we see, in the one-dimensional case, if ( and rj are proportional as functions of ther- 
modynamic quantities such as T and /i, then the bulk and shear viscosity terms are not 
independent, and 

n oc n. 

However, when £ and rj have, in general, different dependence on the thermodynamic quan- 
tities, the two viscosities act differently. 

The time component of the divergence of T^ v gives 

_ e(T) + ^_ + __ = . (50) 

The equation for the space component is automatically satisfied by the scaling ansatz show- 
ing its consistency. The entropy production rate is calculated to be 

^(au" - oaf) = -~ (u + . (51) 

Since II and Q are negative definite, the entropy production is positive definite. 
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A. Solutions 



When ( = £o, f] = 770 and tr are constant, then we can obtain analytic expression for the 
proper energy density. We obtain 



Co 



tr 



a 1 -k) - El -S 1 £ - 



(52) 



where 



Ei(-x) 



oo e - t 



— (it, 
t 



is the exponential integral, and Eq is a constant which should be determined from the initial 
condition for IT (or Q). For relativistic ideal gas, 



we get 



so that for En — 0, 



P = I- 



& , 4e pr/p | \ II (r) 



T \ 4 / 3 



(53) 



(54) 



where the integral containing the exponential function can still be evaluated analytically. 
The temperature is determined from the energy density as 

e = a SB T A , 

where <jsb is the Stephan-Boltzmann constant. 

On the other hand, a typical estimate from the kinetic theory shows that the shear 



viscosity r\ is proportional to the entropy density s, rj = bs, where b is a constant 



Following Ref. 26|, we choose b = 1.1. Furthermore, we use the relaxation time 

_ 3r] IS _ 3?7 
4p 8p 

Here, it should be noted that our definition of rj is twice of other papers 



24 



24 £ | 



(55) 



23|, 



24 
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,1261,12 



The effect of the bulk viscosity has not been discussed in previous papers. We analogously 
assume that the bulk viscosity has the similar s dependence, ( = as. For a baryon free 
relativistic gas, s is related as the energy density as 

s = Ce 3 /\ 
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so that the equation for the energy density becomes 

d 2 e ( 7 t\ de 1/4 



dr 2 



7 t\ de 1 
3 t r ) dr t r \3 



C- 



r 3/4 



T 



(56) 



were 



a 



The above equation is the same as the equation derived in Ref. [26]. It should be noted 



that this coincidence is due to the specific property of this particular model. In Ref. 
the above equation is obtained under the assumption of no acceleration condition, 

d 



dr 



0. 



which is automatically satisfied in the scaling solution. It should be noted that in the IS, 
this condition is not satisfied in general, and one has to solve the original equation of the IS 
without using the above condition. However, our theory does not require such a condition 
at all to be applied. Eq. (156|) can be solved for the initial condition, 



and 



de 
dr 



T=TQ 



e = e (r ) 



- U + n (r ) + \n (r ) 

t 6 6 



(57) 



Now, we show our numerical results. To compare the previous works, we ignore the bulk 
viscosity. In Fig. 1, we show the energy density e obtained by solving Eq. (1561) as function 
of proper time r. As the initial condition, we set s(tq) = 1 GeV/fm 3 , II (r ) = Q (tq) = at 
the initial proper time To = 0.1 fm/c. The first two lines from the top represents the results 
of the LL. The next two lines shows the results of our theory. The last line is the result 
of ideal hydrodynamics. For the solid lines, we calculated with the viscosity and relaxation 
time which depend on temperature fl57|) . Initially, the effect of viscosity is small because 
of the memory effect, the behavior of our theory is similar to that of ideal hydrodynamics. 
After the time larger than the relaxation time, the memory effect is not effective anymore 
and the behavior is similar to the result of the LL. As we have mentioned, the behavior of 
our theory is the same as the result obtained in Ref. in this case. For the dashed lines, 
we calculated with the constant viscosity and relaxation time, 7] = f](eo) and tr = T^(eo)- 
In this case, the viscosity is constant so that the heat production stays longer and has a 
smaller slope as function of time asymptotically. 
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Sometimes the emergence of the initial heat-up in the LL (the dashed curve in Fig.l) 
is interpreted as an intrinsic problem of the first order theory. However, such behavior 
can also appear even in the second order theory. In Fig. 2, we set II (r ) = C( T o)/ T o and 
Q (tq) = i](t )/t as the initial conditions. In particular, the initial heat-up also appear in 
the second order depending on the initial condition for the irreversible currents (see Fig. 2). 
Therefore, this heat-up is not the problem of the first order theory but rather the specific 
property of the scaling ansatz. This was already pointed out by Muronga. The physical 
reason for this heat-up is due to the use of the Bjorken solution for the velocity field. In 
this case, the system acts as if an external force is applied to keep the velocity field as a 
given function of r. Thus, depending on the relative intensity of the viscous terms compared 
to the pressure, the external work converted to the local heat production can overcome the 
temperature decrease due to the expansion. 



VI. SUMMARY AND CONCLUDING REMARKS 



One of the most important questions in the relativistic heavy ion physics is to determine 
the effect of dissipation in flow dynamics of plasma of quarks and gluons. However, a 
consistent calculation of relativistic dissipative hydrodynamics is not trivial at all. Some 
important questions such as propagation of shock wave in the hot QCD matter require a 
careful treatment of the dissipative processes. The Israel-Stewart theory and the divergence 
type theory contain many parameters difficult to be determined from QCD point of view 
and also it is technically difficult to be implemented on practical calculations in its full form 



26 



27|,l39|,Q- 



From the practical point of view, it is not desirable to deal with many of such parameters, 
since the hydrodynamical approach to the relativistic heavy ion reaction processes itself has 
already many uncertainties [4]. In this sense, the IS is not easy to apply at the present stage 
of the analysis of heavy ion physics data. 

Furthermore, when we deal with numerical solutions, we still encounter ambiguities. 
Because dissipations and relaxation times, required to reproduce the bulk properties of fluid, 
are not necessarily the same as those determined from the microscopic theory such as Lattice 
QCD calculations. In a practical numerical calculation, the truncation of hydrodynamic 
degrees of freedom to the discredited variables introduces a natural cut-off frequencies, and 
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the elimination of those degrees of freedom such as small scale vorticities or turbulences will 



appear as effective viscosities 



41] . Furthermore, the standard numerical method for dealing 



with shock phenomena is to introduce the pseudo-viscosity, first invented by von Neumann 



and Richtmyer 



421 ] . It is basically the bulk viscosity depending on the size of the hydro cell. 



For relativistic shock motion, we need to incorporate such mechanism to the hydrodynamical 
code. Considering these aspects, it is desirable to study the effect of dissipation with the 
minimum number of parameters which characterize the physical processes involved and is 
still consistent with the framework of relativity. 

In this paper, we proposed an alternative approach to this question, different from the 
IS. We start from the physical analysis of the irreversible currents according to the Landau- 
Lifshitz theory. Then, the irreversible currents are given by integral expressions which take 
into account the relaxation time. In this way, causality is recovered and at the same time 
a simple physical structure of the LL is preserved. In our approach, only one additional 
parameter was introduced, the relaxation time, tr. The resulting equation of motion then 
becomes hyperbolic and causality can be restored jsj]. Naturally, causality depends on the 
choice of the values of the parameters including the relaxation time. 

perturba- 



35 



36, El 



More specifically, we verified that the linearized equation of motion for smal 
tions in the homogeneous, static background coincides with Hiscock-Lindblom 
except for the coupling among the different irreversible currents. These couplings are not 
included in our theory considering the Curie principle. Of course the Curie principle is 
believed to be valid in the regime of the first order theory and in the second order regime 
these couplings might be present. However, the existence of the Curie principle may imply 
that these couplings are small compared with the direct terms. 

In addition to its simplicity compared to the IS, our formalism has several practical 
advantages. First, the number of independent variables is kept the same as those of the 
ideal hydrodynamics. The irreversible currents are expressed explicitly as memory integrals 
of these independent variables. In extended (irreversible) thermodynamics, the currents are 
created as extended thermodynamic variables and thermodynamic relations are modified 
9), l^j. Thus, in the IS, we should take the modification of thermodynamic relations into 



account but in our case the usual thermodynamical relations should be used. Second, in 
the IS, the projection operators enters in the equations of motion in a complicated way so 
it is not trivial to extract the standard form to apply the conventional method of the time 
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integration. In our case, this does not occur. Instead, we need just some additional integrals 
associated with the local comoving coordinates. Finally, our simplified structure permits 
to include all of the irreversible currents, such as bulk and shear viscosities and thermal 
conduction at the same time without any difficulties. 

The essential difference of our formalism from the IS is the expression for entropy pro- 
duction. In the IS, the entropy production is required to be positive definite algebraically. 
Thus, the integral form like our formulation is not possible even neglecting some coupling 
terms. We relaxed this condition, that is, the positiveness of entropy production is required 
only for the hydrodynamical motion with time scales longer than the relaxation time. For 
extremely violent change of variables, the instantaneous entropy for a hydrodynamic cell 
would not necessarily be positive definite. 

We have applied our theory to the case of the one-dimensional scaling solution of Bjorken 
and obtained the analogous behavior of previous analysis. In this case we can prove explicitly 
the positiveness of entropy production. We showed the time evolution of the temperature. 
As expected, our theory gives the same result of Ref. [2a], because the no-acceleration 
condition used in Ref. 26] is automatically satisfied in this model. Note that our theory is 
applicable to more general case where the acceleration is important. 

The transport coefficients contained in relativistic dissipative hydrodynamics should in 
principle be calculated from QCD. In the first order theory, it is known that the transports 
coefficients can be calculated by the Kubo formula [44]. However, this formula does not 
gives the transport coefficients in the second order theory, as is shown in Ref. 32j explicitly 
and the corresponding corrections should be evaluated. 

Our theory is particularly adequate to be applied to the hydro-code such as SPheRIO 
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22] . Implementation of the present 



which is based on the Lagrangian coordinate system 
theory to the full three-dimensional hydrodynamics is now in progress. 

Authors are grateful for fruitful discussion of C. E. Aguiar, E.S. Fraga, T. Hirano and 
T. Hatsuda. T. Koide would like to thank David Jou for useful comments. This work was 
partially supported by FAPERJ, FAPESP, CNPq and CAPES. 
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FIG. 1: The time evolution of the energy density. The dashed curves correspond to the calculations 
with the constant viscosity and relaxation time. The first two lines from the top represent the 
results of the LL. Next two lines shows the results of our theory. The last line is the result of ideal 
hydrodynamics . 
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FIG. 2: The time evolution of energy density with the different initial conditions from Fig. 1. 
The dashed and short dashed lines represent the result of the LL and our theory, respectively . For 
comparison, our result of Fig. 1 is shown, again (ideal T^^tq)). The last line from the top is the 
result of ideal hydrodynamics. In this case, the energy heat-up is observed even in our theory. 
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